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Abstract 

Einstein's equations are known to lead to the formation of black 
holes and spacetime singularities. This appears to be a manifestation 
of the mathematical phenomenon of finite-time blowup: a formation of 
singularities from regular initial data. We present a simple hyperbolic 
system of two semi-linear equations inspired by the Einstein equations. 
We explore a class of solutions to this system which are analogous to 
static black-hole models. These solutions exhibit a black-hole struc- 
ture with a finite-time blowup on a characteristic line mimicking the 
null inner horizon of spinning or charged black holes. We conjecture 
that this behavior — namely black-hole formation with blow-up on 
a characteristic line — is a generic feature of our semi-linear system. 
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Our simple system may provide insight into the formation of null sin- 
gularities inside spinning or charged black holes in the full system of 
Einstein equations. 
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1 Introduction 

This paper examines a simple system of two equations inspired by the Ein- 
stein equations. The main purpose is to gain insight into the onset of null 
singularities inside spinning or charged black holes (BHs). 

To understand the background and motivation for our toy model it 
will be worthwhile to review the development of our present conception 
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of the null singularity inside BHS. The r = curvature singularity of the 
Schwarzschild geometry has been regarded for many years as a prototype 
for the spacetime singularity expected to be present inside BHs. However, 
the Reissner-Nordstrom (RN) solution, describing a spherically symmetric 
charged BH, lacks a spacelike r = singularity. Instead it admits an in- 
ner horizon (IH) — a perfectly smooth null hypersurface which constitutes a 
Cauchy horizon (CH) for partial Cauchy surfaces outside the BH. The (an- 
alytically extended) RN solution admits an r = singularity too, but this 
singularity is timelike rather than spacelike, and it is located beyond the IH 
(hence outside the Cauchy development). A similar situation is found in the 
Kerr solution, describing a stationary spinning BH: A perfectly smooth IH, 
which again functions as a CH; and the spacetime singularity is timelike, lo- 
cated beyond this null hypersurface. In both the RN and Kerr solutions, the 
regular IH is known to be unstable to small perturbations, and this instabil- 
ity leads to the formation of a curvature singularity instead of a smooth IH. 
Thus, in order to explore the structure of the singularities inside realistic 
spinning BHs, one must understand the process of singularity formation due 
to the instability of the IH. 

Of the three BH solutions mentioned above — Schwarzschild, RN, and 
Kerr — the one which is mostly relevant to realistic spinning BHs is obviously 
the Kerr solution. Nevertheless there is a remarkable similarity between the 
internal structures of spinning and charged BHs, which allows one to use 
spherical charged BHs as a useful toy model for the more realistic (but 
much more complicated) spinning BHs. 

The IH of the RN solution is the locus of infinite blue-shift, as was 
already pointed out by Penrose pp. Infalling perturbations of various kinds 
are infinitely blue-shifted there, which leads to instability of the IH [2]. As 
a consequence the latter becomes the locus of a curvature singularity, to 
which we shall often refer as the IH singularity. In order to explore this 
phenomenon, Hiscock |3] modeled the blue-shifted perturbations by a null 
fluid — a stream of massless particles. He analyzed the geometry inside a 
charged BH perturbed by a single such stream, an ingoing null fluid, using 
the charged Vaidya solution [4]. He found that the IH becomes a non- 
scalar null curvature singularity. Later Poisson and Israel [5] explored the 
system of a charged BH perturbed by two fluxes, namely both ingoing and 
outgoing null fluids. They concluded that in this case too the IH becomes 
a null curvature singularity. This time, however, the singularity is a scalar- 
curvature one because the mass-function — a scalar quadratic in derivatives 
of the area coordinate — diverges, a phenomenon known as mass-inflation. 
The detailed structure of this mass-inflation singularity was later analyzed 
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[6] within a simplified model (in which the outgoing flux is replaced by 
a discrete null shell). This study showed that the metric tensor (when 
expressed in appropriate coordinates) has a continuous and non-singular 
limit at the singularity. Yet derivatives of the metric functions diverge at 
the IH, yielding a curvature singularity. The continuity of the metric has 
crucial physical consequences: It implies that the singularity is weak [7J, 
namely an extended object will only experience a finite (and possibly very 
small) tidal deformation on approaching the IH singularity. 

Subsequently more detailed numerical and analytical studies of the mass 
inflation phenomenon were performed, in which the perturbations were mod- 
eled by null fluids or by a self-gravitating scalar field [3, El QUI CEH 02] • These 
studies confirmed the conclusions of the earlier analyses ([3], [6]). In addi- 
tion, numerical analyses revealed that, at least in the case of scalar field 
perturbations, a spacelike singularity forms in the asymptotically-late ad- 
vanced time. More recently Dafermos |13j proved for a characteristic initial 
value problem for the spherically symmetric Einstein-Maxwell-Scalar Field 
equations that for an open set of initial data on the event horizon (EH), 
the future boundary of the maximal domain of development becomes a null 
surface along which the curvature blows up. Dafermos proved that the met- 
ric can be continuously extended beyond the IH, namely, the singularity is 
weak. 

The situation inside a spinning BH is similar in many aspects to that of 
a spherical charged BH. Here, again, the inner horizon is the locus of un- 
bounded blue shift, suggesting that the regular IH of the Kerr geometry will 
become a curvature singularity when perturbed. A thorough perturbation 
analysis |14j [15] showed that indeed a scalar-curvature singularity forms at 
the early portion of the IH, which is again null and weak. This picture of the 
spinning IH singularity was later confirmed by an independent perturbative 
analysis by Brady et al. [16J. The existence of a class of solutions to the 
vacuum Einstein equations which admit null, weak, scalar-curvature singu- 
larities was also demonstrated in exact non-perturbative analyses \17\ [T8] 
(though these exact analyses, unlike the earlier perturbative analyses, did 
not demonstrate the actual occurrence of a null weak singularity inside BHs). 

If a cosmological constant A > is present, the spacetime is no longer 
asymptotically flat and a cosmological horizon replaces the future null in- 
finity. The spherical charged BH and the stationary spinning BH are then 
described by the Reissner-Nordstrom-de Sitter (RNDS) and Kerr-de Sitter 
solutions, respectively. In both cases there are three horizons, namely cos- 
mological, event, and inner horizons. The surface gravity of these horizons 
depend on the parameters of the solutions, namely the cosmological con- 
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stant, mass, and the charge or angular momentum. We denote these surface 
gravities by k co (cosmological horizon), n ev (EH) and Ki n (IH). If Ki n > k co , 
there is an infinite blue-shift at the IH, suggesting an instability of the lat- 
ter. This instability was first investigated by Mellor and Moss |19| in the 
case of spherical charged BHs, and by Chambers and moss [2U] for spinning 
BHs, using linear perturbations in both cases. In the case of a spherical 
charged BH the non-linear instability with respect to ingoing null fluid was 
investigated by Brady and Poisson [21]. Brady, Nunez, and Sinha [22] inves- 
tigated a model in which both ingoing and outgoing null fluids are present, 
and found that the mass function diverges provided that Ki n > 2k co . In 
the range 2k co > Ki n > k co the mass function is finite, yet the Kretchman 
curvature scalar R a p^sR a ^ 1& diverges at the inner horizon. Later, Brady, 
Moss, and Myers |23j considered also the contribution of the radiation that is 
scattered by the curvature in the vicinity of the EH. They found that when 
this scattering is taken into account, the necessary condition for stability 
of the IH (namely bounded curvature) is that both k co and greater 
than Ki n . None of the stationary (asymptotically-de Sitter) electro- vacuum 
black holes satisfy this condition. Chambers [24] studied a simplified mass- 
inflation model with a continuous ingoing null fluid and a discrete outgoing 
shell and confirmed the earlier results [22] . He also found that for all values 
of Ki n the metric functions are continuous and non-singular at the IH, even 
though the mass function diverges. Namely, the mass- inflation singularity 
is weak in the A > case as well. 

The combination of all the above-mentioned investigations strongly sug- 
gests (though a mathematical proof is still lacking) that the vacuum (or 
electro-vacuum) Einstein equations admit a generic class of solutions in 
which a null weak singularity forms inside a spinning (or charged) black 
hole. In what follows we shall assume that this is indeed the case. Now, 
the Einstein equation in 3+1 dimensions (and with the lack of any sym- 
metry) is a rather complicated non-linear dynamical system. The following 
question therefore naturally arises: Is it possible to extract from the Ein- 
stein equations a smaller and simpler dynamical system, which is capable of 
producing black hole-like configurations with generic null weak singularities 
inside them? If such a simpler system is found, perhaps it could be viewed 
as the "active ingredient" of the Einstein equations (as far as the formation 
of black holes and null singularities is concerned) . This may provide insight 
into the mathematical process of the formation of null singularities. The 
construction and exploration of such a simple system of equations is our 
main goal in this paper. 

The system of Einstein equations combines both evolution and constraint 
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equations. It appears likely, though, that the property of producing generic 
null singularities is admitted by the sub-system of evolution equations. We 
shall therefore extract our simplified toy-system from the evolution equa- 
tions and simply ignore the constraint equations. Now, when the constraint 
equations are discarded, one obtains dynamical behavior even in spherically- 
symmetric situations (it is the constraint section which "freezes" the dynam- 
ics in spherical symmetry). Consequently, we shall extract our toy-system 
from the evolution section of the electro- vacuum Einstein equations in spher- 
ical symmetry. 

We shall thus proceed in Sec. 2 as follows: We start from the electro- 
vacuum Einstein equations in spherical symmetry. We also add a cosmolog- 
ical constant, for reasons explained below. Then we discard the constraint 
equations, and re- formulate the evolution equations in a simple form free 
of first-order derivatives. This yields a semi-linear hyperbolic system of 
two equations for the two unknowns which we denote R(u,v) and S(u,v), 
where u, v are two null coordinate and R, S are constructed from the metric 
functions (specifically R is the square of the area coordinate, and s will be 
specified below). The new system involves a "generating function" h(R) 
which in the above construction emerges in a very specific form [see Eq. 
(|14p ]. However, from the mathematical point of view it appears likely that 
the global properties of the solutions such as BH and singularity formation 
will not be sensitive to the detailed functional form of h(R), but only to 
certain global and/or asymptotic features of this function. For this reason 
we extend our view point and explore this semi-linear system with a rather 
general function h(R). This generalizes our investigation and simplifies it 
at the same time. 

Our strategy of considering a general function h(R) also has a side ben- 
efit: As it turns out, certain two-dimensional general-relativistic dilatonic 
models can be re-formulated such that their evolution sector is described 
by our semi-linear system, with a certain function h{R). This includes the 
model by Callan et al. [25] and its charged generalization [261 G2 EH] . We 
describe this at the end of Sec. 2. 

In Sec. 3 we describe some basic mathematical properties of our semi- 
linear system, including conserved fluxes, a generalized mass function, and 
gauge freedom. The latter means that the semi-linear system is invariant 
under coordinate transformations of the form u — > u'(u),v — > v'(v). Then 
in Sec. 4 we construct, for any h(R), a class of exact solutions with van- 
ishing fluxes, to which we shall refer as the flux-free solutions. This is a 
one-parameter family of solutions (for given h(R)), a generalization of the 
RNDS solution to arbitrary h(R). We then observe that for functions h(R) 
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admitting three roots (or more), the corresponding flux- free solution de- 
scribes a RNDS-like styatic BH, with three horizons, namely three null lines 
of constant R: an event horizon located at a line u = const, and cosmologi- 
cal and inner horizons, both located at v = const. (In the "Eddington-like" 
coordinates, in which the flux-free solution is first derived, the three horizons 
are located at infinite value of the relevant null coordinate, but this is later 
fixed by a coordinate transformation as described below.) The three hori- 
zons intersect at a single point P, representing the timelike infinity for the 
external region between the cosmological and event horizons. The function 
s diverges at P, but this divergence does not represent a spacetime singular- 
ity: Instead it reflects the fact that the proper-time interval between P and 
any point to its past is infinite. The function R is many-valued at P. 

The singularity structure of the flux- free solution is studied in Sec. 5. 
The solution becomes singular at the horizons (s diverges), but this is merely 
a coordinate singularity. To regularize the solution we transform u, v to new, 
"Kruskal-like" , coordinates. (Specifically we "Kruskalize" u with respect to 
the event horizon and v with respect to the cosmological horizon, such that 
the initial data for the BH formation are regular.) In these new coordinates 
the solution extends smoothly into the BH, and provides a description of 
the internal geometry up to the inner horizon. 

The asymptotic form of the functions R, s near the IH is the primary 
objective of this paper. In the flux-free solution (expressed in Kruskal coor- 
dinates and extended into the BH as described above), R admits a constant 
finite value along the IH but s diverges there (for a generic h(R)). This 
divergence, too, does not indicate a true singularity, because it can be re- 
moved by "Kruskalizing" v with respect to the inner horizon. With such a 
coordinate transformation, the variables R, s become perfectly regular (in 
fact analytic) in the IH neighborhood. In fact, this divergence of s reflects 
the infinite blue-shift (or red-shift in some cases) which takes place at the 
IH, just as in the standard RN and RNDS geometries. It should be noted 
that all invariant quantities involving the variables R, s and their derivatives 
are regular at the IH in the flux-free solution. 

Consider now the initial-value problem for our semi-linear system. The 
initial hypersurface is taken to be a spacelike hypersurface which intersects 
both the event and the cosmological horizons. In the first stage we assume 
that the initial data agree with those of the flux- free solution. Then the 

1 Note, however, that such a new "Kruskalization" will spoil the original "Kruskaliza- 
tion" of v at the cosmological horizon, which will be expressed by a divergence of s along 
the latter. 



7 



evolving solution will be just the flux-free solution, with an IH of the form 
described above. We assume that the initial data are everywhere regular, 
which means that the flux-free solution is obtained not in the Eddington- 
like coordinates, but in other coordinates which are regular at the event and 
cosmological horizons — e.g. the above mentioned Kruskal-like coordinates. 

The major challenge is now to understand how will the functions R, s be 
affected if the initial data are modified such that they no longer agree with 
those corresponding to the flux-free solution. What features of the BH and 
the IH will survive the perturbation, and which features will be modified? 
We do not have a full answer to this question, but we do have a conjecture 
that we present in Sec. 6, based on several compelling indications. These 
include the linear perturbation of the flux-free solution, the "generalized 
Vaidya solution" (valid for arbitrary h(R); see the appendix), and also some 
specific examples of h(R) for which the general solution may be constructed. 
Our conjecture may be stated very briefly as follows: First, the global black- 
hole structure is unchanged (this is manifested by the persistent divergence 
of s at a point P where the three horizons meet); Second, R remains finite 
(though no longer constant) along the IH; Third, the divergence of s on 
approaching the IH persists and preserves its leading asymptotic form; and 
after "re-Kruskalization" s becomes finite along the IH, just as in the unper- 
turbed flux-free solution. However, one important difference occurs due to 
the deviation from the flux-free initial data: Although the variables R, s are 
continuous (after "re-Kruskalization") at the IH, they are no longer smooth. 
Certain invariant quantities involving the derivatives of R now diverge on 
the IH (this holds provided that the "surface gravity" of the IH is sufficienly 
large) . 

The h{R) function corresponding to the electro- vacuum solutions in four 
dimensions without a cosmological constant, namely Eq. (|14p with A = 0, 
has only two roots. The corresponding BH solution has two horizons, the 
event and inner horizons. The cosmological horizon disappears when A van- 
ishes, and instead there is a future null infinity. Our semi- linear system is 
useful in this case too, but the initial- value problem described by this system 
is conceptually more complicated in this case. To understand the reason, 
consider a black-hole solution with a cosmological horizon, and consider an 
initial spacelike hypersurface £ which intersects both the event and cosmo- 
logical horizons. We can pick a compact portion So of S which still intersects 
the event and cosmological horizons. Then the early portion of the IH is 
included in the closure of D + (T i q), where D + denotes the future domain of 
dependence. On the other hand, in the analogous asymptotically-flat case 
the initial hypersurface £ must extend to spacelike infinity (or alternatively 
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to future null infinity) in order to have any portion of the IH being included 
in the closure of -D+(S). This means that the behavior of R, s near the IH 
will depend on the asymptotic behavior of the initial data as the initial hy- 
persurface approaches spacelike infinity. No such complication occurs in the 
case of a BH with a cosmological horizon: Here it is sufficient to require that 
the initial data are sufficiently regular on So , and the issue of their large- R 
asymptotic behavior does not arise. For this reason, we shall restrict our 
attention in this paper to functions h(R) with at least three roots. Note that 
in this case the divergence of s at P and at the IH is, from the PDE point of 
view, a manifestation of the finite-time blow-up phenomenon, caused by the 
non-linearity of the hyperbolic system. (The standard General-Relativistic 
point of view is somewhat different, however, because the proper-time dis- 
tance of P is infinite, due to the divergence of s, so this divergence is not a 
spacetime singularity.) 

As was mentioned above, we view our semi-linear system as a toy model 
for the much more complicated system of Einstein equations in four dimen- 
sions. Obviously not all properties of the Einstein equations are mimicked 
by our semi-linear system. The properties we expect our toy system to dis- 
play are (i) the very formation of the BH (expressed in our system by the 
finite-time blow-up of s at the point P), (ii) the no-hair properties of the 
BH — namely the decay of external perturbations, and (iii) the generic for- 
mation of a null, scalar-curvature, weak singularity on the IH. We do not 
expect our toy system to properly address the spacelike singularity (which 
may intersect the IH at later times). Also this simple system is incapable 
of describing the oscillatory character of the null IH singularity inside a 
generically-perturbed spinning BH [29]. Nevertheless our toy system cor- 
rectly demonstrates the basic properties of the IH singularity — in particular 
its weakness. 

We conclude in Sec. 7 by outlining some directions for future research. 

2 The Field equations 

2.1 Maxwell-Einstein equations in spherical symmetry 

We start by considering the Maxwell-Einstein equations in a spherically 
symmetric spacetime. We write the metric in double-null coordinates as 

ds 2 = -2f(u,v)dudv + r 2 (u,v)dn 2 , (1) 
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where d£l 2 = dd 2 + sin 2 9d<p 2 . The Maxwell equations are easily solved, 
yielding 

with all other components vanishing. Here Q is a free parameter, to be 
interpreted as the charge. The electromagnetic energy-momentum tensor 

T,u, = ^ (g^F^F^ - \g^F aP F^ 

is then substituted in the Einstein equations with a cosmological constant 
A, 

Gp, + hgp, = SnT^. (2) 
This yields a system of two evolution equations, 

r, u r, v f ( Q 2 2 \ , , 

r,uv = — 1 - — - Ar (3) 

r 2r \ r z J 

f,uf,v , n / / 2 f 1 r,Q 

and two constraint equations, El 



= ^ + ^, u r, v + ^ 1 - 2^ ) (4) 



^*,mi — T ,uf ,u/ f j — T,vf,v/ f • (5) 

The latter two equations (unlike the two evolution equations) are in fact 
ordinary equations along the lines v = const or u = const, respectively. 



2.2 Constructing our semi-linear system 

We first re-formulate the field equations such that no first-order derivatives 
appear in the evolution equations. To this end we introduce two new vari- 
ables R and s instead of r and /: 

R = r 2 , e s = rf. (6) 

With the new variables, the evolution equations take the convenient form 

Ruv = e s (% + AR? - -L) , (7) 

\R2 Ri) 

2 This involves a slight abuse of the standard terminology as the notion of constraint 
equations is usually formulated with respect to foliations of spacetime by spacelike hyper- 
surfaces. The terminology we use here is a natural extension of the standard one to the 
double-null set-up in two effective dimensions. 
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3Q2 + (B) 



2i?5 2i?2 2R2 
The constraint equations become 

R,uu — R,uS,u j R,vv — R,vS : v • (9) 

In the next stage we simply omit the constraint equations Q and keep 
the evolution equations (|7l8p as our dynamical system. Recall that the evo- 
lution equations form a closed hyperbolic system, which uniquely determines 
the evolution of (properly-formulated) initial data. By this we achieve sev- 
eral goals: First, a non-constrained dynamical system is conceptually simpler 
to analyze than a constrained one; Second, this allows us to explore and test 
our hypothesis that the phenomenon of generic null-singularity formation 
inside four-dimensional spinning BHs is essentially a property of the evo- 
lution sector of the Einstein equations. H Finally, omitting the constraint 
equations retain the dynamics to the problem (it is the constraint sector 
which is responsible to properties like e.g. the Birkhoff theorem). This pro- 
vides us with an effectively two-dimensional toy system aimed at mimicking 
dynamical properties of the Einstein equations in four dimensions. 

Next we recognize that the expression in the parentheses in Eq. (|8j) is 
the derivative of the expression in the parentheses in Eq. ([7]) with respect 
to R. We can therefore write the two equations as 

R, uv = e s F{R) ■ s tUV = e s F'(R), (10) 

where at this stage F(R) denotes the specific function 

F(R) = % + AR 1 2 -^ T , (11) 

R2 R2 

and hereafter a prime denotes a derivative of a function of one variable with 
respect to this variable. These equations, which are semi-linear nonhomo- 
geneous wave equations, constitute the core of our model. 

For later convenience we introduce another function, h(R), defined by 
its derivative: 

F(R) = -h\R). (12) 



3 This is obviously a vague statement because the division of the Einstein equations into 
evolution and constraint subsystems is not unique, but depends on the choice of slicing. 
We expect, however, that this property of the Einstein equations will not be sensitive to 
the details of the foliation chosen. Note that in two effective dimensions the double-null 
formulation induces a unique division into evolution and constraint subsystems. 
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Note that h(R) is defined up to an integration constant. The semi-linear 
system now reads 

R,uv = -e s h'(R) ; s>uv = -e s h"{R), (13) 

In the specific case (fTT|) we have 

i 20^ 2A. 3 

h(R) = 2R^ + —\ i?i + const . (14) 

R2 3 

Note that this is 2r times the standard RNDS function 

1 2m Q 2 2 A 2 

where here the integration constant was expressed as —4m, m being the 
ADM mass of the corresponding static RNDS solution. 

The final stage in constructing our toy system is to abandon the specific 
function h(R) of Eq. (|14p and instead to explore the semi-linear system (|13j) 
for a general function h(R). Again we define F = —dh/dR, or 

h(R) = - j F(R)dR. (16) 

This generalization is advantageous for several reasons. First, if indeed 
the system (|13p leads to generic null singularities, it is plausible that this 
property will not be sensitive to the specific functional form (|14p . Rather, 
we expect the qualitative properties of our dynamical system to depend only 
on certain qualitative features of h(R). Note also that our primary goal is to 
provide a simple toy model aimed at mimicking certain dynamical features 
of e.g. the vacuum Einstein equations in four dimensions, and from this 
perspective the spherically-symmetric electro- vacuum system of the previous 
subsection should itself be regarded as a toy model; hence there is no reason 
to firmly stick to the specific function (I14h . Second, this extension of our 
view-point will allow us to seek simple examples of functions h(R) for which 
the general solution of the system (fT3l) may be constructed. Such solvable 
examples would provide valuable insight into the dynamical properties of 
this system. 



2.3 Application to two-dimensional black holes 

In addition to its role as a toy model for singularity formation, our semi- 
linear system is also directly applicable to certain dilatonic models of two- 
dimensional BHs. In the model developed by Callan et al [25] there is a 
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dilaton (f>(u,v), a cosmological constant of the two-dimensional model A, 
and the metric is 

ds 1 = -e 2p{u ' v) dudv. (17) 

The classical, matter-free, Einstein equations then yield two evolution equa- 
tions and two constraint equations. Transforming to the new variables 
R = e~ 2 ^ and S = 2(p — <fi), the constraint equations reduce to Eq. Q, and 
the evolution equations take the form (I13jl . this time with the generating 
function 

h(R) = X 2 (R + const). (18) 

This dilatonic two-dimensional model was later generalized to include 
a Maxwell field as well as charged matter fields [26\ [271 128j . Here, again, 
with the same substitution R = e~ 2 ^, S = 2(p — <fi) the classical matter-free 
Einstein equations are reduced to Eqs. (|9ll3p with the generating function 

h(R) = X 2 (R + Q 2 /R + const), (19) 

where Q is a parameter proportional to the Maxwell field's charge. 



3 Basic mathematical properties 

In this section we introduce some basic features of our semi-linear system 



3.1 The gauge freedom 

Our semi-linear system (|10p is invariant under a family of gauge transforma- 
tions. These are coordinate transformations which preserve the double- null 
form of the metric: u — * u(u),v — > v(v). The variable R is invariant under 
this coordinate transformation, but s changes. Since e s oc g uv , it transforms 
like a covariant tensor of rank two, and one finds: 

■-=--(£)-Ks)- (201 

The various quantities made of R and s may be classified according to 
the way they transform under a gauge transformation. The simplest are the 
scalars, namely quantities which are unchanged. Obviously R is a scalar. 
Apart from R itself, there is only one scalar made of R and s and their first- 
order derivatives: e~ s R )U R jV . Another useful, non-scalar, quantity is e~ s R tW , 
where hereafter w stands for either u or v. This quantity is invariant to a 
transformation of w, but not to transformation of the other null coordinate. 
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3.2 The conserved fluxes 

Consider the quantities 

3? = R,vv — R,vS,v ; * = R,uu — R,uS,u • (21) 
Differentiation $ with respect to u, one observes that 

\R,uv) a: R,uvS,v R,v8,uv (22) 

identically vanishes by virtue of the field equations (|1U|) . In a similar manner 
one finds that the derivative of *$>(u) with respect to v vanishes. Namely, 

* u = ; *,« = . (23) 

We shall refer to <J?(t>) and ^f(u) as the two conserved fluxes (or simply 
fluxes). It is sometimes useful to express these fluxes as 

$(v) = e s (e~ s R tV ) tV , (24) 
= e s (e-^, u ), u . (25) 

One can easily verify that in a gauge transformation u — * u(u),v —* v(v) 
the two fluxes transform as 

$ = <$>(dv/dvy 2 (26) 
§ = ^(du/du)~ 2 (27) 

(namely like components of a covariant second-rank tensor). 

Note that $(f) and \&(u) are uniquely determined by the initial data 
for R and s (this is most easily seen when the characteristic initial-value 
formulation is used [30]). 

Application to the spherically- symmetric charged case: 
In the four-dimensional spherically-symmetric case, an important prob- 
lem is that of the RN solution perturbed by two fluxes of null fluids, namely 
ingoing and outgoing fluxes. In this case, the dust contribution to the 
energy-momentum tensor is [5j 

T% st = L in (v)/(Anr 2 ) ; 1%? = L out (u) / (Anr 2 ) , (28) 

where r is the area coordinate, u and v are two null coordinates, and Lj n , L out 
denote the ingoing and outgoing dust fluxes, respectively. Prom the Einstein 
equations for T vv and T uu with the line element ([1]) one finds that 

Lin = r {^J^ ~ r ,vv^J ; L out = r f ^' u ^' u - r , uv \ . (29) 
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These quantities are directly related to the conserved fluxes &(v),*$>(u) dis- 
cussed above. In fact one can easily show, using Eqs. (|6|2ip . that 

Lin = -\<$>{v) ; Lout = • (30) 

We may therefore regard the quantities as the generalization of the 

spherically-symmetric null- fluid fluxes to arbitrary h(R). 

It is important to recall that in the spherically symmetric case (jlip the 
semi-linear hyperbolic system (I10p is mathematically equivalent to the mass- 
inflation model [5] with two arbitrary fluxes Lj n and L ou t. 

3.3 The generalized mass function 

Consider the scalar quantity 

M(u,v)=e- s R u R v + h (R) (31) 

where ho(R) is a certain member of the one-parameter family (I16p (namely 
one associated with a certain choice of the integration constant). This is a 
generalization of the mass parameter to dynamical cases (see below). Note 
that for a given F(R) the mass function is defined up to an additive constant 
(associated with different choices of ho (R) ) . 

One can easily show, using the field equation (JT3J) for R, that the deriva- 
tives of M satisfy 

M yV = e~ s R tU <fr(v) ; M u = e~ s R >v ^(u) . (32) 

Also, differentiating the last equation with respect to v and recalling Eq. 
(|24p . one observes that M satisfies the simple field equation 

M >uv = e- s ^{u)<5>{v) . (33) 

From Eq. (I32p we see that when the fluxes $>(v) and vanish, the 

mass function becomes a fixed parameter. This is the situation in the "flux- 
free solution" described below (section [4]) . Also when one flux vanishes (e.g. 
\I / ), the mass function only depends on one null coordinate (v in this case), 
which is the situation in the generalized Vaidya solution (discussed in the 
Appendix). 

Application to the spherically- symmetric charged case: 
As an illustration we consider here the mass function for spherically 
symmetric charged black holes. In this case one naturally defines ho by 
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omitting the constant in Eq. (|14p . namely 



h (R) = 2fiU^-- . (34) 

In terms of the original variables r, / the mass function then reads 

M(u, v) = 2r 1 + ^ - -r 2 + . 35 

In the case of A = 0, the last expression is just four times the function 
m(u,v) defined by Pcisson and Israel 0.0 Rewriting Eqs. O and d33j) 
in terms of r,f,m,a and the fluxes Li n ,L ou t, one recovers Eqs. (4.4) and 
(3.15) therein for m jt) , m jU , and m ini) . 

4 The flux-free solution 

In this section we investigate a class of solutions which is the generalization 
of the static RNDS family to general h(R). These are the solutions in which 
both ^ and $ vanish. We first construct these solutions in Eddington-like 
coordinates and then transform to Kruskal-like coordinates. We then explore 
the singularities of these solutions. 

4.1 Construction in Eddington-like coordinates 

In the case considered here, 

* = $ = , (36) 

Eqs. (I24T25D read 

(e- s i^),„ = ; (e- s R u ) )U = . (37) 
The first integral of these two equations is 

R, v = c u (u)e s ; R )U = c v {v)e s . (38) 

The equation for R v is invariant to a transformation of v, but a transforma- 
tion of u affects the function c u (u). However, the signs of c u is preserved in 
such a gauge transformation, because we require the new null coordinate u 

4 In Ref. [5] e 2<T is used instead of /. 
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to be future-directed, just like the original u. The situation with the equa- 
tion for i? n is exactly the same (with the obvious interchange of u and v). 
Thus, with the aid of a gauge transformation we can bring both functions 
c u (u) and c v (v) to ±1, with the signs corresponding to those of the original 
functions. 

Consider first the case where R is increasing with v and decreasing with 
u, namely c v > and c u < (this is typically the situation outside a BH, 
though no further than the cosmological horizon - namely region I in Fig. 
[I]). Then the gauge transformation described above leads to 




Figure 1: Regions 1,11 and III in the flux- free solution. 

R,v = e s = -R >u . (39) 
This implies that both R and s are functions of a single variable, 

x = v — u , (40) 

and these two functions satisfy 

R, x = e s . (41) 
Then the field equation ([13]) for R uv reads 



R, xx = R, x h'(R) = ^h{R{x)) , (42) 
yielding the first order ODE 

R, x = h(R). (43) 



17 



Thus, R(x) is given by its inverse function 
Then s(x) is given by 

s = ln{h{R)) . (45) 

The second field equation, namely Eq. (p~3|) for s :UV , is automatically satis- 
fied, as one can easily verify. Note that this solution is static, in the sense 
that it only depends on the spatial variable x = v — u. Note also from Eqs. 
(|43|) or (|45|) that h must be positive in this case. 

Next let us consider the case where R is decreasing with both v and u, 
which is typically the situation inside a BH (region II in Fig. [1]). Then 
instead of Eq. ([39]) we now get 

R v = - e s = R yU . (46) 
Correspondingly we now define 

x = v + u , (47) 

and both R and s are functions of x only, satisfying 

R, x = -e s (48) 

this time. Substituting again in the field equations (fl~3[) . one finds that Eqs. 
(|43|) and Eq. (|44|) for x(R) still hold, but there is a sign change in the 
expression for s, namely s = ln(—h(R)). Note that h is negative in this 
case. The general expression for s which holds in both cases is obviously 

s = ln(\h\) . (49) 

Thus, the flux-free solution comes in two versions: The "external-type" 
version, which depends on the spatial variable v — u, and the "internal- 
type" version, which depends on the temporal variable v + u. The first 
version occurs in regions where h > 0, and the second occurs when h < 
(typically inside a BH). Equations (|44|) and ([49 1) hold in both cases. 

5 The full integral of Eq. (|42[) obviously involves two integration constants. The first 
one is already embodied in the definition of h in Eq. (|43p . The second one is an arbitrary 
constant to be added to the right-hand side of Eq. (|44p . This constant has no physical 
meaning, however, as it may be absorbed by a gauge transformation e.g. v —* v + const, 
which merely shifts i by a constant. 
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In its both versions, the flux-free solution is in fact a one-parameter fam- 
ily of solutions (for given F(R)), due to the arbitrary integration constant 
in the definition of h in Eq. (|16p . The situation is somewhat confusing 
because to a specific flux-free solution a function h{R) is associated in two 
different ways: (i) through the field equation (fTHj) (to be satisfied by the 
solution), and (ii) through its explicit construction via Eq. (|44p . It must 
be noticed that h(R) appears in two conceptually different ways in these 
two occasions: In the field equation it appears as an equivalence class (be- 
cause the field equations only depend on dh/dR), and in the construction 
procedure it appears as a one-parameter family of distinct functions. To 
clarify this notational confusion we reserve h(R) of Eq. (I16|) to denote the 
equivalence class, and introduce the notation 

H(R) = h (R) + c , (50) 

where ho(R) is a representative of the equivalence class h(R) (already men- 
tioned above) , and c is an arbitrary constant which distinguishes the various 
members of this equivalence class. Thus, we rewrite the flux-free solution as 

dn 

s = ln(\H(R(x))\) ; — = H(R) . (51) 
(with x = v ± u as above), namely 

x M=fim dR ' (52) 

The parameter c may be interpreted as one associated with the system's 
mass. Indeed, the mass function (j31f) reads for flux-free solutions 

Af (it, v) = h {R) - H(R) = -c . (53) 

In the spherically-symmetric four-dimensional case, the flux-free solution is 
just the RNDS family of solutions, whose ADM mass is related c through 
m = — c/4. We therefore write H in this case as 

i 20^ 2A 3 

H(R) = 2Rk + — x i?i - 4m . (54) 

Ri 3 

4.2 Horizons 

The solutions constructed above become pathological at any value R = Ro 
for which H vanishes. From Eq. (|51|) s diverges there to — oo. Also Eq. 
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([52]) implies (assuming finite F(Rq)) that x diverges at R = Rq, meaning 
that either u or v is unbounded there. This phenomenon is analogous to the 
coordinate singularity at the horizon of the Schwarzschild solution, when 
the metric is expressed in double-null Eddington coordinates. In our case, 
too, this coordinate singularity may be overcome by transforming to new, 
Kruskal-like, null coordinates, as shown in the next subsection. To this 
end, however, we must first analyze the asymptotic behavior of s and x on 
approaching the horizon in the original Eddington-like gauge. 
We define (for each horizon): 

K = H'(R ) (55) 

and assume K ^ 0, therefore 

H(R) ~ K(R - R ) (56) 

in the neighborhood of the horizon. Note that for K > 0, dx/dR = 1/H 
is negative for R < Rq and positive for R > Rq, hence x — > — oo at both 
sides of the horizon. Similarly, for K < at both sides x — > +oo. Therefore 
Kx — > —oo for both K > and K < 0, and at both sides of the horizon. 
In the horizon's neighborhood Eq. (|5ip reads 

R jX ~K(R-R ), (57) 

yielding [f| 

R-R ~ ±e Kx . (58) 
Also the same equation for s reads 

s = ln{\H\) ~ Kx + ln(\K\) . (59) 

Thus, both s and x diverge logarithmically in R — i?o- 

The occurrence of sign flips at the horizon in some of the above expres- 
sions complicates the analysis. To help clarifying this confusion we define, 
for each horizon, the quantity 

X(R) = sign{R - R )e Kx . 

X (unlike x) is continuous and monotonous across the horizon, and it van- 
ishes at the horizon itself. The comparison of Eqs. (j58[) and the definition 

6 In the general solution to Eq. (|57p the right-hand side of Eq. (|58|) should be multiplied 
by an arbitrary constant. However, this constant may be omitted with no loss of generality 
because it can be absorbed by a shift in x. Such a shift merely corresponds to a gauge 
transformation v — > v + const as noted above. 
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of X immediately yields dR/dX = ±1 at the horizon, and a closer look at 
the signs reveals that at both sides 

Since dX/dx = KX, R(X) satisfies the same differential equation at 
both sides: 

dR H(R) 

dx = lcx- (61) 

This, combined with Eq. (|60l) . implies that X(R) is analytic across the 
horizon (provided that H{R) itself is analytic in a neighborhood of R = Rq, 
which we assume). Note that H/X is analytic too, and it gets the non- 
vanishing value 

| = K (62) 

at the horizon. 

We shall primarily be interested in functions H{R) admitting (at least) 
three simple roots Ri (i = 1,2,3), ordered R% < i?2 < R\, such that H 
is positive at R2 < R < R\ and negative at i?3 < R < R2, as shown in 
Fig. [2J An archetype is the function h{R) corresponding to the spherically- 
symmetric electro- vacuum solutions, Eq. (I54p . which (for sufficiently small Q 
and A, and sufficiently large m) admits three roots. These roots correspond 
to the cosmological horizon (i?i), the event horizon (R2), and the inner 
horizon (R3), as shown in Fig. [TJ 




Figure 2: The function H(R), which has three roots: R\, R2 and R3. 

The horizons divide the spacetime into three regions which we denote 
I, II, III, as shown in Fig. [TJ We shall primarily be concerned here with 
the regions I and II (region III will not concern us here, except at the very 
neighborhood of the cosmological horizon). Note that H is positive in region 
I and negative in regions II and III. 

As was shown above, in the Eddington gauge s diverges on the three 
horizons. The divergence at the IH does not pose any difficulty — in fact in- 
vestigating this divergence and its physical implications is one of our primary 
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goals. However, the divergence of s on the event and cosmological horizons 
does pose an undesired feature: We would like to explore a situation in 
which a (flux-free) solution of the type described above, which includes the 
three horizons, emerges from regular initial data prescribed on some com- 
pact spacelike initial hypersurface En. Furthermore we want the horizons' 
intersection point to be included in D + (T,q). To this end, So must intersect 
both the event and cosmological horizons (but not the inner horizon). The 
divergence of s (and also u or v) on these two horizons renders the Eddington 
gauge inappropriate for such a regular initial- value set-up. We shall there- 
fore proceed now to transform the Eddington coordinates into Kruskal-like 
coordinates with respect to the event and cosmological horizons. 

4.3 Transforming to Kruskal-like coordinates 

The construction of the Kruskal-like coordinates in our case is similar to 
the standard procedure in e.g. the Schwarzschild spacetime — except that 
here we "Kruskalize" u with respect to the EH and v with respect to the 
cosmological horizon. 

Let us define on each horizon R = Ri: 

k i =\H'(R i )\ (i = l,..,3) 

(namely, it is the \K\ value associated with the i'th horizon.) Consider first 
the EH, R = i?2- Here K > 0, hence x diverges to — oo. In both regions I and 
II the Eddington coordinate v is regular along the EH but u diverges (see 
Fig. d|). In region I H > 0, hence x = v — u, and the divergence of x means 
that u — > +oo. On the other hand in region II H < 0, hence x = v + u, and 
the divergence of x now implies that u — > — oo. Correspondingly we define 

u = _ e -k 2 u (g3) 

in region I, and 

U = e k2U (64) 

in region II. Then U is continuous across the EH, and is monotonously- 
increasing (namely future-directed) everywhere; it is negative at region I, 
positive at region II, and vanishes at the EH. This transformation cures the 
divergence of s, as we show below. 

7 In the RNDS case, namely the function H(R) of Eq. (|54|l . fci,2,3 correspond to twice 
the quantities n co ^ v ^ n mentioned in the Introduction. 
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Next we consider the cosmological horizon, R = R\. Here K < 0, hence 
x diverges to +00, meaning that v — > +00 on approaching the horizon from 
region I, whereas u is regular. H We thus define (in region I) 

V = -e~ klV . (65) 

Again, V is a future-directed null coordinate which takes negative values in 
region I and vanishes at the cosmological horizon. Since V is continuous 
across the EH, it takes negative values in region II as well. 

The variable R is invariant under the coordinate transformation (u — > 
U,v — ► V). Therefore, R is formally given as a function of U and V through 

R(U,V) = R{x(U,V)) , (66) 

where x(U, V) = v(V) =L u(U), and R(x) is defined through its inverse func- 
tion (|52p . On the other hand s is modified in the gauge transformation 
according to Eq. ([20]) . We shall denote our new Kruskal s by S. It satisfies 

a »du dv 

e=e dujv- (67) 

In both regions I and II we have dV/dv = — k\V, dU/du = ±k2U, and 
e s = ±H. The signs properly combine to yield 

s _ H(R) 

e ~hhuv (68) 



in both regions. Equations (|66p and (|68p constitute the flux-free solution in 
the Kruskal-like gauge. 

We now proceed to show the regularity of R and S, and their smoothness 
as functions of U and V, at both the event and cosmological horizons. Con- 
sidering the EH first, we denote by X2 the function X(R) (defined above) 
associated with the EH. Noting that at the EH K > and hence &2 = K, 
we write 

X 2 = sign(R - R 2 )e k2X . (69) 
One finds (treating carefully the flipping signs) that 

X 2 = -U(-V)~ k21 , (70) 



The detailed behavior of the various functions in Region III will not really concern 
us here. The construction below guarantees that in the Kruskal gauge all the relevant 
functions are analytic across the cosmological horizon, and it will be sufficient for us to 
define these functions through an analytic extension from region I to region III. 
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where hereafter kij = ki/kj > for any i,j = 1...3. Since V is strictly 
negative along the EH, and R is an analytic function of X 2 (as establish 
above for a general horizon), we conclude that R(U, V) is analytic in the 
neighborhood of the EH. To analyze the variable S, we note that 

uv = X 2 (-V) 1+k21 , 

therefore 

• s = 5E^- y )- (1 * > - (71 » 

As was established in the previous subsection, H/X 2 is a regular function of 
R (or X 2 ) which takes the value K = k 2 > at the EH. Therefore S(U, V) 
too is a regular (in fact analytic) function of U and V. 

The cosmological horizon is treated in an analogous manner, except that 
here we do not need to explicitly analyze the various functions in region III: 
Instead we simply extend the relevant functions analytically from region I 
into region III. Note that here K < and therefor k\ = —K. Correspond- 
ingly we find (for region I) 

X 1 = sign{R - R 1 ) e - kix = - e kl{u ~ v) . (72) 

This yields 

X l = V(-U)- kl2 . (73) 

Again we see that R(U, V) is analytic in the neighborhood of the cosmologi- 
cal horizon, because R{X\) is analytic (and U is strictly negative). Recalling 
that 

UV = -Ai(-[/) 1+fcl2 , 

we obtain 

e s = _L Z^g) (_t7)-(i+fci2) . (74) 
k\k 2 X 1 

Again, —hjX\ is a regular function of R (or X\) which takes the value 
— K = ki > at the cosmological horizon. Therefore we conclude again 
that both R(U, V) and S(U, V) are analytic functions in the neighborhood 
of the cosmological horizon. Then in region III R(U, V) and S(U, V) are 
defined to be the analytic extension of the corresponding functions from 
region I across the cosmological horizon. 

In fact it is straightforward to show that R(U, V) and S(U, V) are regu- 
lar at U, V < not only in the neighborhood of the event and cosmological 
horizons, but also in the entire range R 2 > R > R\ and its neighbor- 
hood (provided that H(R) itself is regular throughout). Furthermore, for 
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any smooth spacelike initial hypersurface Eo which intersects the event and 
cosmological horizons, the initial data for R and S (corresponding to the 
flux-free solution in the U, V coordinates) are regular throughout. □ 

5 Singularities in the flux-free solution 

In this section we analyze the singularities that appear in the flux-free so- 
lution using the expressions that were derived in the previous section. As 
before, we shall assume that H{R) has three simple roots at R\, R2 and R3, 
with H positive at R2 < R < R\ and negative at R3 < R < R2 (additional 
roots at R > R\ and/or R < R% are allowed; see last subsection). We divide 
the discussion into two types of singularities — the vertex singularity and 
the IH singularity. 

The vertex singularity 

In the Kruskal-type coordinates defined in the previous section (U, V) = 
(0,0) is the intersection point of the three horizons (see Fig. [1]). At this 
point R is many-valued (for example it admits a fixed value R = Ri along 
the i'th horizon). Furthermore, 

S — > +00 

as U, V — > 0. This divergence proceeds in a slightly different manner along 
various paths towards the vertex (0,0). For example, along curves of con- 
stant R ^ Ri, S oc log \ U\ + log \ V\, as can be seen from Eq. (|68j) . Also, 
Eqs. (HO) and (jUJ) imply that 

e S = J_ (_y)-(i+fc 2 i) ( 75) 
fcl 

along the event horizon, and 

e s = - 1 - (-U)- (1+kl2) (76) 
1^2 

along the cosmological horizon. Since ki2 > 0, S — > +00 along both horizons 
on approaching U = V = 0. 

9 As we show below the solution typically develops a singularity at R = R3, and the 
same applies to roots Rq > Ri of H(R) if such additional roots exist. Nevertheless the 
spacelike hypersurface Eo cannot intersect R = R3 or R = Rq. 
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The singularity at the inner horizon 

Next we analyze the asymptotic behavior of the flux-free solution on ap- 
proaching the IH from region II, namely the limit V — > 0_ in the range 
U > 0. To this end we use the same Kruskal-like coordinates (U, V) intro- 
duced above, which allow for regular flux- free initial data on T,q. 

From Eq. (fBTj) applied to the event horizon we deduce that R decreases 
throughout region II, hence R > R 3 . We denote by X 3 the function X(R) 
(defined in the previous section) associated with the inner horizon. Noting 
that at the IH K < and hence k% = —K, we write for region II (recalling 
X = v + u) 

X 3 = e~ k3X = (-V) kai U~ ki2 . (77) 
Substitution in Eq. ()68[) (which is valid in region II as well) yields 



e s = [_L W. u-Q.+k**)] (_y)^-i. (78) 

kik 2 a 3 



From Eq. (|62p we obtain along the IH H(R)/X 3 = —k 3 . Therefore, for any 
curve in region II which approaches the IH (at some U > 0), the term in 
squared brackets approaches a finite value as V — > 0. Along such a curve e s 
behaves as 

oS ~ ( TT-(i+k 3 2)\ r_v\k 3 i-l 



\hk 2 J 

We find that S — > 00 if k\ > k$, S — > — 00 if k\ < k% and S is regular in the 
case of k% = £3. 

The singularity at the IH is locally gauge-removable; Namely, it may be 
removed by transforming from V to a new Kruskal-type coordinate 

V 3 = -e- k3V = -(-V) k31 , 

defined with respect to the IH. One then obtains in the (U, V3) coordinates, 
in full analogy with Eqs. ([71]) and ([76]) . 



e s 3 = 1 z^^i+ta) (8Q) 
k 2 k 3 A 3 

which admits the regular limit 

e Sz = J_ [/ -(l+fe 3 2) ( 81 ) 

at the IH. However, this transformation will destroy the regularity of S on 
the cosmological horizon. In particular it will spoil the regularity of the 
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initial data on So as the latter intersects the cosmological horizon. Hence 
the divergence of S at the IH (for k\ ^ k%) is globally non-removable. From 
the perspective of the finite-time blow-up phenomenon which concerns us 
here, it is mandatory to introduce regular initial data, and the divergence 
of S on the IH is inevitable. 

We now introduce three additional quantities and explore their asymp- 
totic behavior at the IH. These quantities will serve as gauge-invariant in- 
dicators for certain aspects of regularity or irregularity at the IH. 

i) e~ s Ry R,u '• This quantity is interesting because (unlike S) it is 
a scalar, namely invariant under a gauge transformation. We can therefore 
calculate it directly from the original Eddington-like gauge through Eq. (|51j) : 

e- s R y R,u = e- s {R ;X f = \H\ 

which actually vanishes on the IH. 

(ii) e~ s Ry '• This quantity is not a scalar, yet it is invariant under 
a transformation of V. In a gauge transformation u — > u,v — » v it is 
simply multiplied by du/du. We restrict attention here to transformations 
u — > u which are regular in the interior of region II (i.e. at U > 0). We 
then conclude that although the actual value of e~ s R )V is modified, the 
divergence of this quantity at the IH is invariant to a gauge transformation. 
In the Eddington gauge Eq. (|5ip implies e~ s R jV = — 1. In the Kruskal-like 
gauge we obtain 

e~ s Ry = -dU/du = -k 2 U . 

(iii) S u : This quantity too is invariant under a transformation v — ► v, 
because the term —ln(dv/dv) in Eq. (|20p is independent of U. Therefore 
its regularity at the IH is gauge invariant. We calculate this quantity in the 
Eddington gauge: 

s jU = H >u /H = ^Lr x x !U /H = ^ = -F(R) , 

which is presumably regular. 

The regularity of S t u (or, more generally, of its gauge-transformed coun- 
terpart s^u) is a useful indicator for whether the divergence of S at V = is 
locally gauge-removable or not. To show this, we choose in region II a line 
U = Uq > which intersects the IH. We then define the new coordinate v 
to be the affine parameter along this line, setting e.g. v = at the IH. Then 
s = const = s~q everywhere along this line, and in particular at v — > 0. Now, 
at any U > we have 

-u 



s(U, v) = s + / s, u du . 

JU 
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Therefore regularity of s« implies regularity of s at the IH, and vice versa. 

E 

On the other hand, the divergence of e~ s RyR^u and e~ s Ry indicates 
(qualitatively speaking) to what extent the gauge-invariant variable R is 
smooth at the IH. Recall that after the divergence of S has been removed 
(locally) by a gauge transformation V — > V3, the solution (R,Ss) is con- 
tinuous at the IH, so we are left with the issue of the next-level regularity, 
namely C 1 smoothness at the IH. The quantities e _s RyRjj and e Ry 
both provide information on this issue, in a gauge-invariant manner. 

In the flux-free solution all the above indicators (i-iii) are regular (even 
though S itself diverges). As we shall discuss in the next section, this sit- 
uation presumably changes for indicators (i) and (ii) in the more general 
flux-carrying solutions. 

The discussion above was based on the Kruskal-like coordinates U, V. 
It should be emphasized, though, that the singularity structure is in fact 
gauge-invariant, in the following sense: For any new coordinates U(U), V(V) 
related to U, V in a regular manner, the structure of both the vertex and 
the IH singularities will be the same as in the original coordinates. In other 
words, any choice of null coordinates U, V for which the flux-free initial 
data are regular on So and its immediate neighborhood, will lead to the 
same structure of singularity. (Though, obviously the vertex singularity will 
shift from (0,0) to some other point Uq,Vo.) 

The relation between the PDE and general-relativistic notions 
of singularity 

The issue which primarily concerns us in this paper is that of finite-time 
blowup in a non-linear hyperbolic system of PDEs mimicking the Einstein 
equations. To this end we throughout adopt the standard terminology of 

10 Obviously in the flux- free solutions the local gauge removal of the divergence of S can 
be demonstrated explicitly, see Eq. (|81|l . Yet the indicator S,u is useful for the wider 
class of flux-carrying solutions. 

11 We consider here both quantities (i,h), even though both essentially probe the smooth- 
ness of R at the IH, because each of them has its own advantages and disadvantages. The 
quantity (i) is advantageous because it is a scalar. On the other hand in certain cases 
e~ s Ry diverges, indicating the non-smoothness of R, and yet e~ s R : vR,u vanishes. This 
happens because (qualitatively speaking) although the V-derivative of R diverges, its U- 
derivative vanishes. This situation typically happens in Vaidya-like solutions where only 
ingoing flux is present. In this sense the indicator e~ s R,v is more robust. 

12 by this we also mean that U(U) is monotonic with non-vanishing derivative, and the 
same for V'(V). 
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PDEs in considering the onset of irregularities ("singularities") in the evolv- 
ing solutions. This terminology is quite different from that used in General 
relativity (GR) for discussing spacetime singularities. In this subsection we 
discuss how the various irregularity phenomena admitted by our PDE toy 
system are interpreted from the GR point of view. 

We have seen before that the three-roots flux-free solutions typically 
develop two types of singular phenomena: the vertex singularity and the 
IH singularity. It should thus be emphasized that both phenomena do not 
imply a General-relativistic spacetime singularity, as we now discuss. 

The divergence of S at the vertex U = V = fails to mark a General- 
Relativistic singularity because (in a typical application of our system to a 
general-relativistic spacetime) no divergence of curvature is involved. (Note 
that S is a gauge-dependent quantity, and gauge-invariant quantities like 
e~ s R yR u do not diverge at the vertex.) The divergence of the metric 
function g uv cx e s simply implies that the proper-time distances between 
the vertex and points in region I are infinite — namely, the vertex becomes 
future timelike infinity. 

Similarly, the divergence of S at the IH fails to mark a spacetime singu- 
larity because it may be locally removed by a gauge transformation V — > V3. 
Smoothness indicators like e _s RyRu or e ~'~ 'R,v (and others not mentioned 
here) all indicate that subsequent to such a coordinate transformation the 
manifold is perfectly regular (in fact analytic) at the IH. 

Despite the failure of the vertex and IH singularities to mark genuine 
GR singularities, both phenomena actually bear crucial implications to the 
structure and features of the GR black-hole spacetime, as we now discuss. 

The vertex singularity functions as timelike infinity, as we already men- 
tioned. In fact it defines the event and cosmological horizons — both are null 
lines which intersect the vertex singularity at their future. El Furthermore, 
from the standard hyperbolic-PDE point of view the occurrence of a singu- 
larity at the vertex implies the presence of a Cauchy horizon in the (U, V) 
manifold at V = 0. This property carries over to the GR terminology — the 
inner horizon at V = is indeed a Cauchy horizon of the spacetime. Thus, 
the very presence of a Cauchy horizon in e.g. the RNDS family of black 
holes may be attributed to the divergence of S at the vertex. 

Next let us consider the General-relativistic implications of the IH sin- 
gularity. As was mentioned above the divergence of S at V = may be 

13 In the flux-free solutions the horizons may alternatively be defined as null lines of 
constant R, or as the locus of H — 0. These definitions do not hold in the more general 
flux-carrying solutions. 
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gauge-removed in a local manner only: Curing the divergent S at the IH 
will inevitably lead to blow-up of S at the cosmological horizon. The mis- 
match of S between the inner and cosmological horizons is an inherent, 
gauge-invariant, global feature of the typical three-root flux-free solutions. 
In the general-relativistic language this is translated into infinite blue-shift 
which takes place at the IH of the RNDS black holes (as well as their ro- 
tating counterparts). This divergent blue-shift, in turn, leads to a null 
curvature singularity at the IH when the BH is generically perturbed. But 
this relativistic inner-horizon singularity is represented in our PDE termi- 
nology by the divergence of quantities like e~° RyRjj and e~ s ' Ry (which, 
as was mentioned above, only diverge when fluxes are added). 

Three roots versus four roots 

The analysis and discussion so far was independent of the behavior of H(R) 
at R > R\. The extensions of the spacetime diagram in Figs 1 or 3 to the 
right of the cosmological horizon will depend the properties of H(R > R\). 
We may conceive of several options: 

(i) H(R) admits an additional root at Rq > R%; 

(ii) No additional root exists at R > R\; either (iia) H(R) diverges 
(or admits an irregularity) at a certain finite R s > R\\ or (iib) H{R) is 
regular — and negative — for any finite R > R\ (the asymptotic behavior of 
this function as R — > oo will not concern us here). 

Consider first the case (i). In this case the dynamics in region III is 
essentially the same as in region II. In particular, the line in Fig. 3 
(the IH) which borders region II has a counterpart at the future boundary 
of region III — namely the null line located at U = 0, V > 0, which we 
denote Hq. The behavior of R and S on approaching Hq just reflects that 
at H3, with the obvious interchange of the horizons' indices 3 — > and 
1 — ► 2,2 — > 1. In particular R gets the fixed value Rq, whereas S diverges. 
The sign of divergence is —00 (implying infinite blue-shift) for ko > &2 and 
+00 (infinite red-shift) for ko < &2- Essentially the spacetime diagram of 
Fig. 3 (with the Hq singularity added) would exhibit a reflection symmetry 
with respect to a vertical line which passes through the vertex at U = V = 0. 
The singularity will have a symmetric "V-shape" in this case, namely two 
null singularities which emerge from a common vertex singularity. 

Note that additional roots at R > Rq will not affect the spacetime di- 
agram because we terminate the solution at the -Ho singularity, namely at 

14 This holds in the case £3 > fci. In the other case kz < fci there is an infinite red-shift 
at the IH. 
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R = i?o- For the same reason, additional roots at R < R$ will not affect the 
spacetime diagram of Fig. 3 (regardless of whether a root Rq > Ri exists 
or not). 

Turn now to discuss the case (ii), namely H{R) has no roots at R > 
R\. In this case region III will terminate at a spacelike boundary line, 
corresponding to R = R s in case (ha) or R — > 00 in case (iib). This spacelike 
line, which intersects the vertex U = V = 0, will constitute a singularity 
of our PDE system. Recall, however, that in the GR terminology the 
spacelike boundary will not necessarily be a spacetime singularity; Instead 
it may well mark a regular boundary at timelike infinity. 

The function H(R) of Eq. (|54p corresponds to the RNDS black-hole 
solutions. In this case there are only three roots, corresponding to the cos- 
mological, event, and inner horizons. This function belongs to the class (iib), 
and the future spacelike boundary of region III marks the future timelike 
infinity of the external de Sitter universe. 

6 The conjecture about the singularity formation 

The semi-linear system (|13|) is the main objective of this paper and the sub- 
sequent paper [30] . In the previous sections a class of special solutions was 
constructed which demonstrate a finite time blow-up. We saw various fea- 
tures of the singularity which evolved from regular initial data. The solutions 
that we constructed were flux-free, namely "static" . They all correspond to 
flux-free (\P = $ = 0) initial data. This leaves us with the following open 
question: How will these solutions change if we perturb them by modify- 
ing the initial data, such that (^,$ 7^ 0)? We are particularly interested 
in the finite-time blowup phenomenon and the structure of the singular- 
ity exhibited by the flux-free solutions: To what extent it is stable to the 
perturbations? 

To be more specific, and to set the notation for the discussion below, let 
us consider a specific generating function h(R) and a specific three-root (or 
more) flux-free solution H(R) = ho(R) + c. The flux-free initial data are 
prescribed on an initial hypersurface So which intersects the cosmological 
and event horizons (where R = R\ and R = R2, respectively). The initial 
data need not correspond to the Kruskal-like coordinates U, V, but could 
be in any other coordinates U, V which are regular functions of respectively 
U and V — namely, we only require the flux-free initial data to be regular 
through So- The evolving solution then develops a vertex singularity at 

15 S diverges in case (iia), and R (and possibly also S) diverges in case (iib). 
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P = (Uq,Vq), and a null IH singularity at V = Vq. Now we add small 
(but finite), regular, perturbations to the initial functions on So such that 
7^ 0). Which properties of the flux-free singularity will survive the 
perturbations and which will be modified? 

We conjecture that the basic structure of a vertex singularity and a null 
IH singularity is rather robust and will survive the perturbation (see Fig. 
3). However, some of the features of the IH singularity will be modified. 




Figure 3: The general singularity structure. The three horizons are displayed, 
denoted by Hi (i — 1, 2, 3). The singular IH is displayed by a thick line emerging 
from the vertex singularity (the intersection point of the three horizons). 

In particular, the character of the flux-free IH singularity of being fully 
locally-removable (by means of a local gauge transformation) will no longer 
hold when fluxes are added. A more detailed description of the robust and 
non-robust properties of the flux-free singularity is given below. 

Our conjecture is based on several sources of evidence: (i) A linear per- 
turbation analysis [31]. Note that at the linear level the two fluxes ^ and 
do not interact. The effect of each flux may be separately explored by 
means of the exact Vaidya-like solutions (see Appendix A), and then these 
two perturbations are simply superposed. This drastically simplifies the 
analysis, as no PDE needs be solved. In addition, the perturbation scheme 
suggests that the effects of non-linear coupling between the two fluxes are 
negligible compared to the basic effect of the blue-shifted incoming flux. 
(This insight — namely the dominance of the linear perturbation at the IH 
over the non-linear effects — has emerged in previous studies of the IH singu- 
larity inside perturbed spherical charged BHs [Bjlllj. and, more generally, of 
four-dimensional general-relativistic null spacetime singularities |14^ \T7\ [18] ; 
and it applies to our toy semi-linear system as well.) (ii) The general exact 
solution of our semi- linear system for the case h = cos(R) |31j . 

This case is exactly solvable, and it demonstrates all the elements of 
the conjecture — except that it satisfies = k\ hence it does not test issues 



32 



related to the blue-shift phenomenon. Nevertheless it provides support to 
our conjecture concerning all properties of the vertex singularity, as well as 
some of the properties of the IH singularity — in particular its null character, 
(hi) Mathematical study of the general, flux-carrying solutions in the case 
of a saw-tooth function h(R) [30J. In this case h(R) is made of three linear 
sections, patched together to form a continuous but non-smooth function. 
This type of h(R) confirms all elements of the conjecture below — for both 
cases &3 = ki and k^^k\. 

In what follows we present in detail our conjecture about the singularity 
structure when the flux-free solution undergoes a generic small perturbation. 
We divide the discussion into the vertex singularity and the IH singularity. 

6.1 The perturbed vertex singularity 

The perturbed solution exhibits a finite-time blow-up and develops a vertex 
singularity at a slightly-modified point P' = (Uq, Vq), 0. where 

• s diverges logarithmically to +00 on approaching P' (e.g. from the 
past); 

• The function R is many-valued. 

We stress again that this divergence of s does not represent a spacetime 
singularity — it merely reflects the fact that the proper time interval between 
P' and any point to its past is infinite. Therefore P' functions as timelike 
infinity in the BH spacetime. The presence of the singular point P' defines 
the horizons, and in this sense the very formation of the black hole. More 
specifically, the cosmological and event horizons are now defied to be the 
two null lines which intersect P' at their future [namely, (U < UL V = Vq) 
and (U = Uq, V < Vq), respectively]. These two null lines are no longer lines 
of constant R. Nevertheless we conjecture that 

• along each of these lines R admits a well-defined value R[ (i = 1...2) 
on approaching P'; 

• Furthermore these two values are both roots of the same function 
H'(R) = ho(R) + d, with d close to the original c (hence R[ 2 are close 
to the two original roots R\,2)- 

16 Throughout this section we use a prime to denote perturbed quantities and not the 
derivatives of functions 
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E 

From the GR point of view, this conjectured structure of the perturbed 
vertex singularity is a manifestation of the no-hair principle to our toy sys- 
tem; Namely all initial perturbations (the non-vanishing fluxes) should dis- 
perse away and leave a black-hole spacetime which at late time is well de- 
scribed by the corresponding static solution. 

6.2 The perturbed IH singularity 

First we conjecture that, when a flux- free solution with an inner horizon is 
weakly-perturbaed, no spacelike singularity will form in the neighborhood 
of P'. Namely, in the neighborhood of P' the perturbed solution will be 
regular throughout V < Vq, U > Uq. However the solution will generically 
develop a null singularity at the IH, namely at V = Vq (for U > Uq). This 
singularity will be similar in many respects to its flux- free counterpart: 

• R has a well-defined finite limit Rjh(U) at V — > Vq; 

• Furthermore as U — > Uq (the P' limit) Rw{U) approaches a limiting 
value -R3 which too is a root of H'(R), close to the original root R3. 

• Just as in the unperturbed case, the function s generically diverges on 
the IH. The sign of divergence is determined by the relation between 
k[ and k' 3 , where k[ is \dH' /dR\ at the root R = R\: Namely s — > — 00 
for k' 3 > k[ and s —* +00 for k' 3 < k[. El These two cases correspond 
to infinite blue shift or infinite red shift, respectively. 

• In the special case k[ = k' 3 we expect s to remain regular at V — ► Vq. 
In this case the IH becomes a perfectly regular Cauchy horizon of our 
nonlinear hyperbolic system. 

• Returning to the generic case k[ ^ k' 3 : Although s diverges, this 
divergence is locally removable by a gauge transformation 

V 3 == -(V^-V) k z /k K 

17 In fact, the same property is shared by R' 3 , namely the limiting value of R on ap- 
proaching P' along the perturbed IH, as discussed in the next subsection. 

18 If the original flux- free solution satisfies £3 < k\ or k-j, > ki, then we may assume 
that the perturbation is sufficiently small — and correspondingly c' is sufficiently close to 
c — such that the same inequality is satisfied by the perturbaed quantities k' 13 . Obviously 
in such a case, in the above condition for the sign of the divergent s we may replace k[ 3 
by ki,3. 
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But this transformation leads to a divergence of s along the cosmo- 
logical horizon. Namely the divergence of s at the IH is globally non- 
removable. 

• The indicator ds/dU is regular at the IH limit V — > Vq. As discussed 
in the previous section, this indicates that indeed the divergence of s 
is locally removable. 

All the above features merely reflected the similarity between the per- 
turbed and unperturbed (flux-free) IH singularities. But there also is an 
important difference between the perturbed and unperturbed singularities 
in the case of infinite blue shift. This difference is manifested by the fol- 
lowing properties which all hold for generically-perturbaed solutions with 
k' 3 > 2k[ l 19 l(to which we may refer as "strongly-divergent blue shift"): 

• The quantity e~ s (dR/dV) diverges at the IH limit V — ► Vq. 

• The same applies to the scalar quantity e~ s {dR/ dV){dR/ dU). 

• As a direct consequence of the above the mass function, which was 
fixed (M=-c) in the flux-free case, diverges at the IH of the generically- 
perturbed solution. 

The meaning of these last indicators is simple: The divergence of s at the IH 
may be locally removed by a gauge transformation as mentioned above. This 
transformation renders the solution continuous at the IH. Yet generically the 
flux-carrying solution fails to be smooth at the IH. This lack of smoothness 
is manifested by the divergence of certain gauge-invariant quantities, so it 
cannot be cured by any gauge transformation. In the case of strongly- 
divergent blue shift (k' 3 > 2k[) the solution fails to be C 1 . As it turns 
out, in the case k[ < k' s < 2k[ the solution is C 1 but fails to be C 2 at 
the IH. In the case k' 3 < k[ (infinite red-shift) the solution will be C 2 , yet 
generically smoothness will fail at a certain derivative orders (unless 2k[/k'^ 
is an integer). [^| 

In both blue-shift cases (k' 3 > 2k[ and k[ < k' s < 2k[) the generally- 
perturbed solution fails to be C 2 at the IH. In General-Relativistic terms 
this means that the Reimann tensor will diverge there (e.g. as measured by 
polynomial scalars or by its p.p. components) — even if the mass function 
remains bounded. On the other hand, in the red-shift case (k' 3 < 2k[) the 

19 This condition was derived in the perturbed RNDS case in Ref. [22] . 
20 The general rule is the following: The conjectured smoothness level of R as a function 
of V 3 is the same as that of the function \V 3 \ 2k 'i/ k 3 . 
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Riemann tensor will remain bounded, though generically its derivatives will 
blow up at a certain order. 

Our primary concern here is the case of generic perturbations, in which 
both ingoing and outgoing fluxes are present. However, it is also interesting 
to consider the case a BH perturbed by an ingoing flux while the outgoing 
flux vanishes. This case is exactly solvable by means of the Vaidya-like solu- 
tion (see Appendix). In this case scalar quantities as M and e~ s (dR/dV)(dR/dU) 
remain bounded. However e~ s (dR/dV) generically diverges at the IH if 
k' 3 > 2k[. Although this quantity is not a scalar, its divergence is a gauge- 
invariant phenomenon as discussed in the previous section. Again this di- 
vergence indicates an unbounded curvature (though only in the p.p. sense). 

As was mentioned above, the divergence of s at the IH is conjectured 
to be removable by means of a local gauge transformation, leading to a 
continuous (though not smooth) solution at the IH. From the GR point 
of view this means that the metric tensor has a continuous (though not 
smooth) non-singular limit at the IH. This means that the IH singularity is 
weak (namely physically non-destructive), as discussed at the Introduction. 

We summarize some of the conjectured properties of the generically- 
perturbed singularity in the table below. 



The quantity 


Type of divergence 


Condition for 










the divergence 




no null fluid 


only influx 


two null fluids 






+00 


+00 


+00 


k\ > &3 


s 








k\ = &3 




— 00 


—00 


—00 


k 3 > ki 


e~ s R, v 




+00 


+00 


k 3 > 2ki 


e R tV R tU 






+00 


k 3 > 2ki 



Table 1: The divergence of various quantities at the inner horizon H3 expected to 
be found in various cases. 



7 Future research directions 

The problem which motivated this paper concerned the internal structure 
of generically-perturbed spinning vacuum black holes. One would like to 
prove (or disprove) the generic formation of a null singularity inside the BH, 
and its weakness. We find it likely that our simple two-dimensional semi- 
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linear system presented here properly mimics these properties of General- 
Relativistic spinning black holes. 

Prom this perspective the task of mathematically exploring the internal 
structure of realistic spinning black holes may be divided into the following 
stages: (i) Proving (or disproving) the conjecture presented in the previ- 
ous section concerning the generic behavior of our semi-linear system; (ii) 
Extending the analysis to the asymptotically-flat case, namely to functions 
H(R) with two roots and with large-i? asymptotic behavior similar to the 
A = case of Eq. (JMJ; and (iii) Further extending the analysis to the real 
problem of a generically-perturbed spinning black hole. 

Although stages (i) and (ii) are not easy, stage (iii) will probably pose a 
much harder challenge, as it involves the transition from two to four dimen- 
sions (and with a much larger number of unknown functions). Nevertheless 
it may be hoped that the insight gained from our simple toy system in stages 
(i,ii) may make this challenge a bit easier. 

Acknowledgements: We thank G. Wolansky for many discussions and 
comments on the manuscript. In addition we thank J. Feinberg for com- 
ments on a preliminary version of the manuscript. 

A The case of one flux (Vaidya-like solution) 

In the case of a flux in one direction only our semi-linear system may be 
reduced to a single ordinary differential equation along null lines. Consider 
for example the case $ = 0, namely outflux only (the other case \P = may 
be treated in a fully analogous manner). Then Eq. (pHJ) yields R tV = c u (u)e s . 
After a re-labeling of u this becomes 



where we have considered here the case of positive c u for concreteness (the 
case c u < proceeds in a similar manner). Substituting this in the right- 
hand side of the first evolution equation in (fT3j) we obtain 



e s = R 



(82) 



R 



,uv — 



R tV h'(R). 



Integration with respect to v gives 



R >u = -h (R) + M u (u) 



(83) 



where M u (u) is an arbitrary function of u. 
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Differentiation of Eq. (|82p with respect to u, substitution of the evolution 
equation for R and then a second differentiation with respect to v, one finds 
that the second evolution equation (fT3j) is satisfied as well. 

Thus, the original system (fT3j) of evolution equations has been reduced to 
a single ordinary differential equation (|83p along lines of constant v. Solving 
this ODE (with a v-dependent initial condition Ro(v) at a certain initial u 
value u = uq) yields the function R(u, v). Subsequently s is obtained by Eq. 

Substituting this solution in Eq. (|3"Tj) one obtains 

M(u,v) =M u (u) . 
The outflux is then given by Eq. (|25|) or (1321). 



(and, recall, $ = 0). 



T , , dM u , ,, 

• („) = (84) 
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